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ASYMPTOTICALLY  ROBUST  DETECTION  OF  STOCHASTIC 
SIGNALS  IN  CONTAMINATED  NOISE 

M.S.  SCHNIT7.ER  and  D.R.  HALVERSON 
Department  of  Electrical  Engineering 
Texas  A  &  M  University 
College  Station,  Texas  77843 

ABSTRACT 

We  consider  the  discrete  time  detection  of  stochastic 
signals  in  white  noise,  where  the  univariate  noise  density 
is  known  perfectly  only  on  an  interval  about  the  origin. 

We  present  a  method  to  enhance  the  asymptotic  performance 
of  the  detector  by  exploiting  this  knowledge,  and  at  the 
same  time  preserve  robustness  properties  of  the  detector 
to  the  remaining  inexact  knowledge  of  the  univariate  noise 
density  via  a  saddlepoint  condition.  We  then  provide 
examples  to  show  that  improved  performance  is  indeed 
obtained. 


_ INTRODUCTION 

Consider  the  discrete  time  detection  of  stochastic 
signals  in  independent  and  identically  distributed  noise. 

It  the  underlying  statistical  distributions  are  completely 
known,  then  selection  of  a  Neyman-Pear son  optimal  detector 
is  possible.  However,  in  many  cases  knowledge  of  the 
univariate  noise  density  is  imperfect,  particularly  for 
extreme  values  of  its  argument.  In  such  situations  the 
employment  of  a  detector  which  is  robust  to  the  inexact 
statistical  knowledge  may  be  desirable.  Earlier  work  of 
Kas sam  and  Thomas  (lj  and  Poor,  Mami,  and  Thomas  [2]  points 
toward  the  introduction  of  limiting  for  observations  of 
large  magnitude  in  order  to  inp art  robustness  into  the 
detection  scheme.  In  particular,  [2]  considers  the 
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asymptotically  robust  detection  of  stochastic  signals  in 
imperfectly  known  noise. 

In  this  paper  we  consider  further  applications  of 
detectors  possessing  the  same  property  of  limiting.  We  note 
that  often  more  is  known  about  the  univariate  noise  density 
near  the  origin  than  on  the  tails,  and  such  knowledge  has 
the  potential  to  be  exploited  for  improved  performance.  In 
some  earlier  work  [3]  we  show  how  to  design  the  robust 
detector  by  exploiting  this  knowledge  for  improved  performance 
over  the  detector  of  [1];  this  work  is  applicable  to  the  case 
of  a  time  varying  deterministic  signal  embedded  in  imperfectly 
known  noise.  We  now  consider  the  analogous  extension  of  [3] 
for  improved  performance  over  the  detector  of  [2].  Our  results 
illustrate  how  such  improved  performance  can  be  obtained  and 
still  preserve  the  robustness  of  the  detector  to  the  remaining 
inexact  knowledge  of  the  univariate  noise  density  through  the 
employment  of  limiting. 

II.  DEVELOPMENT 

Suppose  we  denote  the  i.i.d.  discrete  time  noise  process 
by  {N.;i=l,2,  ...}  ,  with  f(«)  denoting  the  associated 

univariate  noise  density,  and  let  the  stochastic  signal  be 
denoted  by  { S j ; i =  l , 2 , . . . }  ;  the  zero  mean  independent 

signal  process  is  assumed  to  be  independent  of  the  noise.  The 
detection  problem  then  reduces  to  a  choice  between 
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where  we  observe  realizations  y^  of  the  random  variables  Y^. 
As  in  |l-3],  we  assume  the  noise  density  f  (  •  )  is  a  member  of 
the  Huber-Tukey  mixture  class,  i.e.  f  (  •  )  can  be  expressed  as 
f(.)=(l_i)f0(.)+L.h(.),  where  t  is  a  fixed  known  number 


between  zero  and  unity,  fg(«)  is  a  known  symmetric  nominal 


each 


x  on  an 


density,  and  h(*)  is  an  unknown  "contaminating"  density.  In 
addition,  however,  wc  assume  that  f(x)  i (ATSC) 
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interval  [-a,  a|,  i.e.  h ( x )  =  f^(x)  lor  all  x  in  [-a,  a).  We 

note  in  passing  that  a  similar  assumption  is  employed  in  [2] 
without  being  directly  exploited.  Motivated  by  the  results  of 
(1-3J,  as  well  as  the  form  of  the  locally  optimal  detector 
given  by  [2],  we  employ  the  class  of  detectors  where  a  cest 
statistic  of  form 

Tg  (  y  )  =  g(y.)  is  compared  to  a  threshold,  with  the 

i  =  l  1 

nonl inearitv 


g( • )  given  by  g(x)  = 


(x)/fQ  (x) 


f  "  ( k ) / f  (k) 


x  |  <k 


otherwise 


We  thus  see  that  specifying  k  effectively  specifies  the 
detector. 

In  order  to  analyze  and  compare  the  performance  of  the 
detector,  and  hence  determine  k  for  best  performance,  an 
asymptotic  fidelity  criterion  is  employed.  Because  of  its 
relationship  to  the  Pitman-Noether  theorem  [4],  the  criterion 
of  Asymptotic  Relative  Efficiency  (ARE)  is  particularly  useful; 
the  associated  efficacy  functional  n  provides  a  convenient 


measure  of  performance,  where  n  is  defined  by 
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in  other  work  which  employs  the  ARE,  the  results  of  this  paper 
are  valid  only  when  various  mild  regularity  conditions  are 
imposed  on  the  underlying  statistical  distributions.  These 
assumptions  are  imposed  to  guarantee  that  our  expressions 
exist  and,  when  necessary,  are  nonzero;  additionally. 


'""‘erro 


regularity  assumptions  art-  required  to  insure  the  relevancy  of 
the  efficacy  functional  to  the  ARE  via  satisfaction  of  the 
hypothesis  of  the  i’  i  t  ma  n- Noe  t  her  theorem  |  4  ]  .  In  this  paper, 
our  results  are  presented  under  the  tacit  assumption  of  these 
numerous  though  obvious  mild  conditions.  We  refer  the 
interested  reader  to  |5-7],  in  which  analogous  assumptions  are 
explicitly  listed.  It  should  be  noted  that  we  do  not  impose 
symmetry  assumptions  on  the  contaminating  density  h(*)- 

To  preserve  the  robustness  property  of  the  detector,  k  is 

then  chosen  according  to  sup  inf  n(k,h),  where  dependence  of  r\ 

k  h 

on  k  and  h(-)  is  indicated.  Such  a  choice  of  k  will  be  called 
an  optimal  choice.  The  following  result  provides  insight  into 
how  to  obtain  an  optimal  choice  of  k: 


Theorem:  Let 


K  -  {  k>0  :  if  k>a  then  |g  (  •  )|  and  |g  '  (  •  )|  are  weakly  monotoneton 

la,k]}  .  Then  a  necessary  and  sufficient  condition  for  k  to 
be  an  optimal  choice  over  K  is  that  k  yield  the  supremum  among 


(a)  inf  n(k,h)= - " — 
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(b)  inf  p,  (  k  ,  h  )  *  - 

h  k 
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it  k  a  . 


Proof :  Consider  first  the  proof  for  part  (a).  Note  that  by 

negative  scaling  when  necessary  we  can  assume  without  loss  of 


l  .  e  . 


generality  that  g'(.)  is  weakly  increasing  on  |a,k|, 
g"lx)  £0  on  [a,k|.  Moreover,  by  symmetry  we  also  have 
g"(x)>0on[-k,-a].  in  a  manner  similar  to  the  proof  of  the 
Theorem  of  [ 3 ]  we  obtain  by  direct  evaluation 

0  (  k  ,  h  )  = 


g"(x) fQ(x)dx  +(l-i  ) 


g'*(x)f0(x)dx 


+  t 


_  a 


g" (x)h (x)dx 


+  g"(x)h(x)dx 
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(g2  (x)-g2  (k  )  )  f()(x)dx  +2(1-0  f  (g2  (x  )-g2  (k))f0(x)dx 
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+g-(k)  +t 


-a 


(g2(x)-g2(k))h(x)dx  +  |(g2(x)-g2(k))h(x)dx 


Noting  the  weak  mon  o  t  on  ic  i  t  y  To  f  |g(.)|  implies  that  g2  (  •  )  is 
weaklv  increasing,  we  observe  that  the  numerator  of  the  above 
is  minimized  while  the  denominator  is  simultaneously  maximized 
if  h  ( x )  =0  on  [_k,  _a]  l>  la,k].  The  expression  given  for 


inf  n(k,h)  in  part  (a)  then  follows  directly, 
h 


For  part  (b)  a  similar  evaluation  when  k£a  yields 

l  0 


n  (  k  ,  h  ) 
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Noting  that 


h  (  x  )  d  x  + 

h(x)dx=  1-2 

.  . 

• 

(  x  )  dx  , 


we  obtain  an  expression  for  n(k,h)  which  is  independent  of 

h ( • )  and  which  agrees  with  the  expression  given  for  inf  n ( k , h ) 

h 

i  n  pa  r  t  (  b  )  . 


We  remark  that  differentiation  of  the  expressions  of  the 

Theorem  could  be  employed  as  in  [3]  for  the  purpose  of 

obtaining  a  condition  on  an  optimal  choice  k  as  a  solution 

of  a  nonlinear  equation.  However,  we  have  found  that  for 

this  stochastic  signal  situation  the  appropriate  equations 

are  typically  so  complicated  that  it  is  no  more  difficult 

to  search  directly  over  the  inf  q  (  k  ,  h  )  values.  We  also 

h 

remark  that  the  techniques  of  [3]  could  be  employed  to  extend 

the  Theorem  to  the  case  of  an  asymmetric  nominal  density;  the 

chief  cost  for  such  an  extension  would  be  the  considerably 

more  cumbersome  expressions  for  the  inf  q(k,h)  quantity. 

h 

We  now  consider  an  example  which  employs  a  zero  mean 
unit  variance  sech  noise  nominal  density,  i.e. 
f0  ( x )  =  l"2  sech(4rx)  for  all  x. 


Wt  compare  the  "worst 


case"  efficacy 


2  w 


of  the  detector 


of  | 2]  (i.e.  the  minimal  efficacy  over  those  h(*)  satisfying 

our  constraints  as  well  as  those  of  [2])  with  the  "worst  case' 

elficacy  1,  of  our  detector  with  k  chosen  via  the  Theorem. 

1  ,w 

The  results,  which  depend  on  the  choice  of  a  and  t,  are 


tabulated  below: 


T able: 

Values  of  • 

1  w/  '2  w 
*  * 

.=0.1 

,=0.3 

>  =0. 5 

O 

II 

a  =  2 . 00 

6 . 02 

2.65 

1  .  A  4 

1.13 

a- l . 50 

6.0  2 

2 . 65 

1.44 

1.13 

a=l .00 

3.45 

2.7  3 

1.44 

1.13 

O 

tn 

o 

II 

2  .  73 

1.55 

1.15 

1  .  04 

Note  that  the  most  dramatic  improvement  in  performance 
occurs  when  the  greatest  certainty  is  present  concerning  the 
univariate  noise  ciensitv  (small  t  and  large  a).  This  is 
probably  consistent  with  what  we  might  expect,  however  it  is 
also  directly  opposite  Lo  what  was  observed  for  the  Gaussian 
example  of  [3]  for  the  time  varying  deterministic  signal  case. 
Other  such  opposing  effects  are  known  to  exist  when  comparing 
the  stochastic  and  deterministic  signal  cases,  e . g .  respective 
symmetry  and  skew  symmetry  of  the  nonlinearity  g(-)>  and  thus 
the  results  illustrated  here  need  not  be  viewed  as  surprising 
when  compared  to  those  of  [ 3  J . 

M I.  CONCLUSION 

We  have  introduced  an  approach  toward  improving  the 
asymptotic  performance  of  a  detector  of  a  stochastic  signal 
by  exploiting  available  knowledge  of  the  univariate  noise 
density  near  the  origin,  while  at  the  same  time  retaining  the 
robustness  of  the  detector  to  inexact  knowledge  of  the  tails 
ot  the  noise  density.  We  have  shown  by  way  of  example  that 
such  improved  "worst  case"  performance  does  indeed  occur. 
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